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1 Introduction 

In this paper, we will consider the Benjamin-Bona-Mahony (BBM) [Q like equation (B(m,n)) with a fully 
nonlinear dispersive term of the form 

ut + u x + a (u rn ) x - {u n ) xxt = 0, m, n > 1, m ^ n . (1) 

This equation is similar to the nonlinear dispersive equation K(m, n), 

u t + (u m ) x + (u n ) xxx = 0, to > 0, 1< n < 3 (2) 

which has been studied in detail by P. Rosenau and J.M. Hyman (2). In the literature there are many studies 
dealing with the travelling wave solutions of the K(m,n) and B(m,n) equations, but in general they are 
restricted to the case ttt, = tt, E] H SI |5] |6] |7] El |9l [10] [II] [12 [13 [141 . When m ^ n, the solutions of K (m, n) 
were investigated in [2, 3 |. Our aim here is just to search for solutions of the equations B(m, n), with to ^ n, 
by means of the factorization method. 

We remark that this method |[T5l[T6l[T7l[T8l . when it is applicable, allows to get directly and systematically 
a wide set of solutions, compared with other methods used in the BBM equations. For example, the direct 
integral method used by C. Liu f\3j can only be applied to the B(2, 1) equation. However, the factorization 
technique can be used to more equations than the direct integral method and also, in some cases, it gives rise 
to more general solutions than the sine-cosine and the tanh methods ifPTl [T2l l20l . This factorization approach 
to find travelling wave solutions of nonlinear equations has been extended to third order nonlinear ordinary 
differential equations (ODE's) by D-S. Wang and H.LHI2T1. 

When we look for the travelling wave solutions of Eq. (Q3, first we reduce the form of the B(m, n) equation 
to a second order nonlinear ODE and then, we can immediately apply factorization technique. Here, we will 
assume m 7^ n, since the case m = n has already been examined in a previous article following this method 

ma. 

This paper is organized as follows. In section 2 we introduce factorization technique for a special type of 
the second order nonlinear ODE's. Then, we apply straightforwardly the factorization to the related second 
order nonlinear ODE to get travelling wave solutions of B(m, n) equation in section 3. We obtain the solutions 
for these nonlinear ODE's and the B(m, n) equation in terms of Weierstrass functions in section 4. Finally, in 
section 5 we will add some remarks. 
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2 Factorization of nonlinear second order ODE's 



Let us consider, the following nonlinear second order ODE 



(3) 



where (3 is constant and F(W) is an arbitrary function of W. The factorized form of this equation can be 
written as 



T9 -h(w,o) 



W(9) = 0. 



(4) 



Here, /i and /ja are unknown functions that may depend explicitly on W and 9. Expanding and comparing 
with (01, we obtain the following consistency conditions 



_ F(W) d/i 

71/2 w d# 



TO 
/2 + ^w^- /3 - 



If we solve <(5j for j\ or /2, it will supply us to write a compatible first order ODE 



d6 



W{6) = 



(5) 



(6) 



that provides a solution for the nonlinear ODE 01 lfT31 [161 [171 [181 . In the applications of this paper f% and 
will depend only on W. 



3 Factorization of the BBM-like equations 

When Eq. (0 has the travelling wave solutions in the form 

u[x,t) = 4>(£), £ = hx + wt 



(7) 



where h and w are real constants, substituting into (JTJ and after integrating, we get the reduced form of Eq. 
to the second order nonlinear ODE 



AS- Bi 



Notice that the constants in Eq. (0 are 



A = 



h + w 
h 2 w 



B = 



hw ' 



D = 0. 



D = 



R 

h 2 w 



(8) 



(9) 



and R is an integration constant. Now, if we introduce the following natural transformation of the dependent 
variable 

0"(O = W(6), £ = 9 (10) 

Eq. ® becomes 

- - • a.) 



d9 2 



AW~ -BW~ +D = 0. 



Now, we can apply the factorization technique to Eq. dTTb . Comparing Eq. 0) and Eq. dTTb . we have (3 — 
and 

F(W) = -(AW* +BW- -D). (12) 
Then, from 01 we get only one consistency condition 



fl + fiW^-AW 1 -^ - BW^ + DW- 1 = 



(13) 
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whose solutions are 



h{W ) = ±±J^W^ + ^-W^-2DW + C (14) 
W V n + 1 m + n 

where C is an integration constant. Thus, the first order ODE © takes the form 



^J^Aw^ + ^± W ^. 2DW + C ^. <1 5 > 
dx) V + l m + n 

In order to solve this equation for W in a more general way, let us take W in the form W = tp p , p ^ 0, 1, then, 
the first order ODE (15[ is rewritten in terms of ip as 

,dtp. 2 2nA „/i-»< | n 2nB / m-» s , 9 2D 2 __ C 2 _ 2 „ 

^ = 27 ZTT ^ 5 + 27 o_ " ' ~ — V 2 P + -V P - (16) 

ax) p (n + 1J p z (m + n) p z 

If we want to guarantee the integrability of ( fTST ), the powers of <p have to be integer numbers between and 
4 ll24l . Having in mind the conditions on n, m (n ^ m > 1) and p (p ^ 0), we have the following possible 
cases: 

• If C = 0, D = 0, we can choose p and m in the following way 

2n n+1 3n- 1 

p = ±- with m=— — , — - — ,2n-l (17) 

1 — n 2 2 

and 

77- 

p = ±- with m = 2n-l,3n-2. (18) 

1 — n 

It can be checked that the two choices of sign in dTvT > and (fT~8b give rise to the same solutions for Eq. (Q]i. 
Therefore, we will consider only one of them. Then, taking p — — 1^377-. Eq- ( Q becomes 

,dip. 2 A(n-l) 2 B{n-l) 2 n+l 
(S) 2 = ^7^ + «;^ (20) 



d6> ; 2n(n + l) n(5n- 1) r ' 2 

dp 2 = A(n-1) 2 B(n-\) 2 
de' 2n(n + l) n(3n-l) 



,dw, 9 A(n-l) 2 Bfn-l) 2 

(^) 2 = „ \ 1 + m = 2n-l (21) 



and for p — — -r^- . 



,dip, 2 2A(n~l) 2 2B(n-l) 2 „ 
a# n(n+l) n(3n— 1) 

,du> s0 2A(n-l) 2 B(n~l) 2 , 

J = ) mU 77 ^ ) m = 3n-2. (23) 

aw n (n + 1) n(2n — 1) 

If C = 0, we have the special cases, p — ±2, n = 2 with m = 3, 4. 

Due to the same reason in the above case, we will consider only p = 2. Then, Eq. ( fl~6l ) takes the form: 

,dw. 9 DA B n „ 

(J) 2 -- y+ 3^ +y ^ m = 3 (24) 
,d</\ 2 D A B 4 , 
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• If A = C = 0, we have p = ±2 with m = — , > 2 n - 
In this case, for p = 2, Eq. ( TTo*i > has the following form: 



W = ~2^ + 3^ ' m= 2 (26) 

(|)--f^ + f, — 



• If A — 0, we have p = ±1 with m = 2 ri, 3 n. 

Here, also we will take only the case p = 1, then, we will have the equations: 



(^) 2 = -2Dip+ ^-Bip 3 + dp\ m = 2n (29) 
at) 3 

(^) 2 = - 2 ^V+f ^ 4 + C, m = 3n. (30) 



• If A = D = 0, we have p = ±— with m — 3n,5n. 

2 

Thus, for p = — , Eq. ( fT~6b becomes: 

A 2 = 2 Bp 3 + 4CV, m = 3n (31) 
(^) 2 = ^V 4 + 4C^ m = 5n. (32) 

4 Travelling wave solutions for BBM-like equations 

In this section, we will obtain the solutions of the differential equations (fT9b-(l23l in terms of Weierstrass 
function, p(8; 52, 53), which allow us to get the travelling wave solutions of B(m, n) equations (HJ. The rest 
of equations (I24l i-(l32l can be dealt with a similar way, but they will not be worked out here for the sake of 
shortness. 

First, we will give some properties of the p function which will be useful in the following l25ll26ll . 

4.1 Relevant properties of the p function 

Let us consider a differential equation with a quartic polynomial 

(S) 2 = P M = a °^ 4 + 4ai<^ 3 + 6a 2 ip 2 + 4a 3 tp + a A . (33) 
do 

The solution of this equation can be written in terms of the Weierstrass function where the invariants g 2 and 173 
of (O are 

52 = a o a 4 ~ 4 ai 03 + 3 a\ , g$ = a 2 &4 + 2 a\ a 2 0,3 — a 2 — flo a 3 — a \ a 4 (34) 
and the discriminant is given by A = g\ ~ 27 gf. Then, the solution cp can be found as 



ip(6) =tp + jPA'Po) (p(0; 32,53) - j^-FWOpo) ) <35) 
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where the subindex in P^, (<£>□) denotes the derivative with respect to ip, and ipo is one of the roots of the 
polynomial P((p) ( f33l >. Depending of the selected root ipo, we will have a solution with a different behavior 

mm. 

Here, also we want to recall some other properties of the Weierstrass functions |27|: 

i) The case g 2 = 1 and g 3 = is called lemniscatic case 

p(9;g 2 ,0)=g 1 2 /2 p(9g 1 2 /4 ;l,0) 1 g 2 > (36) 

ii) The case 172 = —1 and g 3 = is called pseudo-lemniscatic case 

p(6;g 2 ,0) = \g 2 \ 1 / 2 p(6\g 2 \ 1 / 4 ;-l,0), g 2 <0 (37) 

iii) The case g 2 = and 53 = 1 is called equianharmonic case 

P(9;g 2 ,0) = g 1 3 /3 p(9gl /6 ;0 1 l), g 3 > . (38) 

Once obtained the solution W(9), taking into account (0, dTOb and W = ip p , the solution of Eq. (Q~|i is 
obtanied as 

u(x,t) = = W"{6) = <p%(6), 9 = £ = hx + wt. (39) 



4.2 The case C = 0, D = 0, p = 

n + l 
• m = 

2 

Equation ( fT9l can be expressed as 



2n 
1 — n 



l d6» ; V ^ 2n(n + l) n(3n+l)^ 
and from P(y) = 0, we get the root of this polynomial 

^(3n + l) 

<A) = 7 ■ (41) 

V 2B(n + l) ' 

The invariants (|34l ) are: g 2 — g 3 — 0, and A = 0. Therefore, having in mind p(0; 0, 0) = -Jr, we can 
find the solution of ( fT9] l from ( |35l > for <^o, given by ( HTt . 

P 2 (n-l) 2 (n + l)fl 2 -2An(3n+l) 2 
^ = 4Pn(n+l)(3n+l) ' (42) 

Now, the solution of Eq. (Q]) reads from d39l 

"P 2 (n - l) 2 (n + l) (hx + wt) 2 - 2 An (3n + l) 2 



3n- 1 



4Pn(n + 1) (3n+ 1) 



(43) 
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In this case, our equation to solve is (|20| > and the polynomial has the form 

A(n~l) 2 P(n-l) 2 3 



2 n (n + 1) ?i (5 ?i — 1) 



5 



with one real root: (po 
ants 



-A (5 

2 B (n+1) 



.92 = 0, 



1/3 



Here, the discriminant is different from zero with the invari 

-AB 2 {n-lf 



93 



32 n 3 (n + l)(5n-l) 2 ' 
Then, the solution of ( f20b is obtained by ( 1351 ) for <y9o, 

"4n (5 n - 1) p(0; 0, g 3 ) + 2B (n — l) 2 ^ 



An (5 n - 1) p(0; 0, ,g 3 ) - B (n - l) 2 
and we get the solution of Eq. (QJ from (f39b as 



u(x, i) 



4n(5n - I) p{hx + w t;0, g 3 ) + 2 B (n - I) 2 ip 
4 n (5 n — 1) p(/i x + w t; 0, 33) — B (n — l) 2 (fo 



(45) 



(46) 



(47) 



_ ( -A (5w-i; 
2 B (n+1) 



1/3 



with the conditions: A < 0,g 3 > 0, for ip = ( 

the solution d4Tb in terms of equianharmonic case of the Weierstrass function 



Using the relation fl38l >, we can write 



u(x, t) 



-A(5n-1) \ 2/3 ( 2 2 / 3 p{{hx + wt) gl /6 ;0,l) + 2 
2B(n + l) / \2 2 / 3 p((hx + wt)gl /6 ; 0,1) - 1 



(48) 



m = 2 n — 1 

In Eq. ( TSTT i. the quartic polynomial is 

P(p) = 



A(?i-1) 2 B{n-lf 



and has two real roots: tpo = ± 
invariants are 



-A (3ti-1) 
B (n+1) 



2ra(ra + l) 2n(3n-l) 
1/4 
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(49) 



32 



for A < 0, B > or A > 0, B < 0. In this case, the 
AB(n- l) 4 



.93 = . 



(50) 



4n 2 (n + l)(3n-l)' 

Here, also the discriminant is different from zero, A ^ 0. We obtain the solution of ( EH from d35l l for 

¥>o, 

- 4n {n + 1) ip 2 p{0; g 2 ,0) - A (n- If 



9 ^° [An(n+l)ip 2 p(0;g 2 ,O)+A(n-l) 2 
and we get the solution of Eq. (fTJ from d39b as 



u{x, t) 



2 , An (n + 1) (fl p(hx + wt;g 2 , 0) - A(n — 1) ! 
r " ' An(n+ l)iplp(hx + wt;g 2 ,0) + A(n - If 



(51) 



(52) 



with the conditions for real solutions: A < 0, B > 0, gi < or A > 0, B < 0, 52 < 0. 

Having in mind the relation d37l i. the solution d52l can be expressed in terms of the pseudo-lemniscatic 
case of the Weierstrass function: 



u{x, t) 



-A(3n- 1) 



B(n+ 1) 
for A < 0, B > 0, g 2 < and 

' -A (3ri-l) 



it(a;, £) 



B(n+1) 



1/2 



1/2 



2p((fox + ^)| g2 | 1 / 4 ;-l,0) + l 
2p((/ l x + W i)l<? 2 | 1/4 ;-l,0)-l 



2p((fox + u>Q| g2 | 1 / 4 ;-l,0)-l 
2p((/ l x + u-i)|<? 2 | 1 / 4 ;-l,0) + l 



(53) 



(54) 



for A > 0, B < 0, g 2 < 0. 
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4.3 The case C = 0, D = 0,p 

• m = 2 n — 1 

Now, the polynomial is cubic 



n 



1 — n 



2A{n-l) 2 2B{n-l) 2 



n (n + 1) 



■<P- 



n (3 n — 1) 



(55) 



and has three distinct real roots: ipo = and ipo = ± 
0. Now, the invariants are 



-A (3n-l) 
B(n+1) 



1/2 



for A < 0, B > Oor A > 0, B < 



.92 



-AB (n - l) 4 
n 2 (n+ 1) (3 n- 1) : 



.93 = 



and A^O. The solution of ( f22l is obtained from ( f35l > for i/jq, 

pn(n+l)^ p(^;g 2 ,0)-A(n-l) 2 
^ ^° [2n(n+l)<p Q p(0;g 2 ,O)+A(n-l) 2 

and substituting (l57b in d39l >, we get the solution of Eq. (Q]l as 

'2 77 (n + 1)990 p(^a; + wt;g 2 ,0) -A(n- l) 2 



i) 



2n(n + l)(p p(/ix + wt;32,0) + A(n- 1) : 



with the conditions: A < 0, B > 0, .g 2 > and A > 0, B < 0, 52 > for <^ 



-A (3n-l) 
B(n+1) 

-A(3n-1) 
B (n+1) 



(56) 



(57) 



(58) 

1/2 
1/2 



While the root ip = leads to the trivial solution, u(x, t) = 0, the other root (p = 
gives rise to imaginary solutions. 

Now, we can rewrite the solution d58l in terms of the lemniscatic case of the Weierstrass function using 
the relation ( f36l > in d58l i: 



u(x, t) 



A(3n- 1)V /2 / 2 p((/i a; + wt)g 2 '' 1 ; 1,0) + 1 



1/4.- 



+ / V 2 p(( /ia; + wf )ff2 1/4 ; 1 ) o) - 1, 



(59) 



for A < 0, 5 > 0, g 2 > and 



u(x, t) 



-A{Zn-l) 
B(n + l) 



1/2 



' 2p{{hx + wt)g\ ,A ;l,Q)-l 
2p((hx + wt)g 1 2 /4 ;l,0) + l / 



(60) 



for A > 0, B < 0, g 2 > 0. 

771 = 3 71 — 2 

In this case, we have also a quartic polynomial 



P{}p) 



2 A (77- l) 2 5(77- l) 2 j 



It has two real roots: <£>o = and = ( 

.92 = 0, 



?7(77+l) ' <y3 + 77 (2 77 - 1) 

1/3 



(61) 



2A(2n-l) 
B (n+1) 



.93 



For the equation (l23l l. the invariants are 

A 2 B(77- l) 6 



4n 3 (t7 + 1) 2 (2t7- 1) 



(62) 
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and A^O. Now, the solution of d23l reads from (l35l l for cpo, 

2n(n + l)ip p(d;0,g 3 )-A(n 



<P = <Po 



2n{n+l)<p o p(0;0,g 3 ) + 2A{n-l) 2 
Then, the solution of Eq. (Q~|i is from d39i > as 

2 n (n + 1) ifo p(h x + w t; 0, 33) — A{n— if 



u(x, t) 



¥0 



(63) 



(64) 



k 2 n (n + 1) ip p(h x + w t; 0, y 3 ) + 2 A (n - l) 2 / 

with the conditions: B < 0, 33 > 0. Taking into account the relation (1381 . this solution also can be 
expressed in terms of the equianharmonic case of the Weierstrass function: 



u{x, t) 



2A(2n-l) 
B(n + 1) 



1/3 



1 



2 2 / 3 p((hx + wt)gl /6 ;0,l 
2^p((hx + wt)g 1 3 /e ;0,l) + 2 



(65) 



We have also plotted these solutions for some special values in Figs. (fx)-©. We can appreciate that for the 
considered cases, except the parabolic case ( 142b . they consist in periodic waves, some are singular while others 
are regular. Their amplitude is governed by the non-vanishing constants A, B and their formulas are given in 
terms of the special forms (|36ll-(l3~8ll of the p function. 
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Figure 1 : The left figure corresponds to the solution d54l ) for h 
right one corresponds to the solution d53l for h = l,w= 1, a = 



-2, w — 1, a — — 1, n — 3, m = 5 and the 
-1, n = 3, m = 5. 
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Figure 2: The left figure corresponds to the solution d54b for h 
right one corresponds to the solution d53l for /i = 1, tu = 1, a = 



-2, u; = 1, a — — 1, n — 2, m = 3 and the 
-1, n = 2, m = 3. 



5 Lagrangian and Hamiltonian 

Since Eq. ( ITTb is a motion-type, we can write the corresponding Lagrangian 



At 



2 n - 



B n m +n 

w— 



m + n 



DW 



(66) 
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Figure 3: The left figure corresponds to the solution d59b for h 
right one corresponds to the solution d60l ) for h = 1, w = 1, a - 



—2, w = 1, a = — 1, n = 3, m = 5 and the 
— 1, n = 3, m = 5. 
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Figure 4: The left figure corresponds to the solution (f59b for /i 
right one corresponds to the solution ( f60b for /i = 1, w = 1, a = 



-2, io = 1, a = — 1, n = 2, m = 3 and the 
-1, n = 2, m = 3. 



and, the Hamiltonian = WgPw — L\y reads 



fr„r(w,iv,fl) = 

where the canonical momentum is 



P 2 

1 M 



n + 1 



m + n 



2L>T^ 



(67) 



(68) 



The independent variable 6* does not appear explicitly in d67| i, then H\y is a constant of motion, Hw = E, with 

2 



1 



(dWY 



2 An Tir n+i 2Bn TJ .m+« 



n + 1 



m + n 



2L>W 



(69) 



Note that this equation also leads to the first order ODE dl3T > with the identification C — 2 E. Now, the energy 
E can be expressed as a product of two independent constant of motions 



where 



n + 1 m + n 

and the phase 5(6') is chosen in such a way that I± (6) be constants of motion (dl± (9) /d6 = 0) 

AWi +BW^ - D 



S{9) = 



:d6. 



2 An W S±l + 2Bn W ^ 



(70) 



(71) 



(72) 



n+1 



2L>iy 
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Figure 5: The left figure corresponds to the solution d48l > for h = —2, w = 1, a = — 1, n = 2, m = 5/2 and 
the right one corresponds to the solution do3T > for h = l,w = 1, a = — 1, n = 3/2, m = 5/2. 

6 Conclusions 

In this paper, we have applied the factorization technique to the B(m, n) equations in order to get travelling 
wave solutions. We have considered some representative cases of the B(m, n) equation for m ^ n. By using 
this method, we obtained the travelling wave solutions in a very compact form, where the constants appear 
as modulating the amplitude, in terms of some special forms of the Weierstrass elliptic function: lemniscatic, 
pseudo-lemniscatic and equiaharmonic. Furthermore, these solutions are not only valid for integer m and n 
but also non integer m and n. The case m = n for the B(m, n) equations has been examined by means of 
the factorization technique in a previous paper [ 14 1 where the compactons and kink-like solutions recovering 
all the solutions previously reported have been constructed. Here, for m ^ n, solutions with compact support 
can also be obtained following a similar procedure. We note that, this method is systematic and gives rise 
to a variety of solutions for nonlinear equations. We have also built the Lagrangian and Hamiltonian for the 
second order nonlinear ODE corresponding to the travelling wave reduction of the B(m, n) equation. Since the 
Hamiltonian is a constant of motion, we have expressed the energy as a product of two independent constant 
of motions. Then, we have seen that these factors are related with first order ODE's that allow us to get the 
solutions of the nonlinear second order ODE. Remark that the Lagrangian underlying the nonlinear system 
also permits to get solutions of the system. There are some interesting papers in the literature, where starting 
with the Lagrangian show how to obtain compactons or kink-like travelling wave solutions of some nonlinear 
equations l|28l |29l [30l [3T1 l32l . 
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